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Abstract 

We study the robustness of complex networks subject to edge removal. Several network models 
and removing strategies are simulated. Rather than the existence of the giant component, we use 
total connectedness as the criterion of breakdown. The network topologies are introduced a simple 
traffic dynamics and the total connectedness is interpreted not only in the sense of topology but 
also in the sense of function. We define the topological robustness and the functional robustness, 
investigate their combined effect and compare their relative importance to each other. The results 
of our study provide an alternative view of the overall robustness and highlight efficient ways to 
improve the robustness of the network models. 

PACS numbers: 89.75.Fb, 89.75.Hc 



* Electronic address: 



lisheng@sjtu.edu.cn 



I. INTRODUCTION 



Complex networks are ubiquitous in our world. They exhibit not only diverse structural 
characteristics HQ , such as the power-law tail of degree distribution and the small-world 
phenomenon of average path length, but also different levels of robustness, e.g., scale- free 
networks display higher tolerance to error but more vulnerability to attack than exponential 
networks jf],!?, S, 9|. Robustness evaluates the ability of a network to maintain its original 
attributes and functions when constituent loss or other kinds of damage are present. Once 
the robustness is identified, the weakness of a network is pointed out for optimizations 
or countermeasures, e.g., the weakness of a communication network can be overcome to 
increase reliability, while the weakness of an epidemic network can be utilized for efficient 
destruction. For these practical applications, the study of the robustness have received a lot 
of interests. 



Several measures of the robustness 



existence of the giant component 



rave been proposed. A frequently used one is the 



m 

13]. 



, Is], Q|- The giant component is the only component 
in a network whose size scales linearly with the number of vertices. It was found that the 
damage to a network, such as random removal of vertices or edges, can be exactly mapped 
to a standard percolation process. The network percolates if the giant component exists, 
indicating that the general connectedness of the network is maintained. There are also other 
_ ha S ed on pities such as efficiency Q, Q and avera g e path ,en gt h fl ~ 
The measure we adopt in the present work is the preservation of total connectedness 
The total connectedness of a network is preserved only when communication is effective 
between every pair of vertices in the network. The more damage needed to destroy the 
total connectedness, the more robust the network is. This measure can be useful to describe 
the robustness of networks that have no vertex redundancy. In such networks, each vertex 
contributes to the whole in a way that cannot be replaced by the others and even the 
unavailability of a single vertex affects the overall functionality, e.g., in a scenario that 
many computers collaborate through a network to accomplish some calculation-intensive 
task, if one or more members lose communication, the performance of the collaboration 
may be degraded or the task may even fail. Note that the total connectedness used in the 
measure is not limited within the sense of topology, i.e., physical connection is not the only 
factor that influences the effectiveness of the communication; there are other factors that 
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can hinder the communication even though all the vertices are connected. These factors 
are often related to the function of a network, e.g., it was reported that in a transmission 
network, apath connecting two vertices may become so long after damage that this path is 



unusable 



16j |. As both topology and function related factors can prevent the communication 



between vertices, two types of robustness are involved: the topological robustness and the 
functional robustness. Studying their combined effect and the dominance of one over the 
other can provide additional insights into the network robustness. However, such studies 
have not yet been carried out extensively. 

Besides the measures, the damage of constituent loss has also been widely studied in 
the form of vertex removal. Removing a vertex is an appropriate abstraction of several real- 
world events, e.g., a user leaves a P2P network, a website goes offline. There are other events 
that should be abstracted more appropriately by removing an edge, e.g., a network cable is 
unplugged, a flight between two cities is canceled. However, only a few works cover edge 
removal. In fact, there are remarkable differences between the two types of removal. In the 
sense of topology, vertex removal inflicts more damage, as each time a vertex is removed, all 
of its edges are removed as well. Thus to achieve the same effect, smaller amount of vertices 
is removed than edges. If the damage of removal is assessed in terms of edges rather than 
vertices, the vulnerability of a network could be moderated jl2j ]. In the sense of function, 
the removal of vertices reduces the total amount of quantity transmitting on a network, as 
vertices usually not only deliver but also generate the quantity. If the reduction of the total 
amount can compensate the damage, the robustness of the network is enhanced 13]. On 
the other hand, the enhancement may not be expected for edge removal because there is 
no change in the total amount. These differences indicate that edge removal needs separate 
studies from vertex removal. 

In this paper, we propose to study the effect of edge removal on the topological robustness 
and the functional robustness, using the total connectedness as the criterion of breakdown. 
Several network topologies and removing strategies are simulated. We combine a simple 
traffic dynamics with the network topologies and successively remove edges until the total 
connectedness of the networks is destroyed. The fractions of removed edges characterize the 
robustness of the networks. Moreover, as the destroy of the total connectedness is related 
to either the topological factors or the functional factors, we see which type of the factors is 
dominant. The purpose of our study is to explore the relationship between the topological 



3 



robustness and the functional robustness, which is different from the purposes of previous 



works that also deal with edge removal 
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221 ] . Particularly, we only concern 



the times when the networks lose the total connectedness, and do not discuss the evolutions 
thereafter, such as the occurrence of cascading failure ll| . The results of our study provide 
an alternative view of the network robustness and highlight efficient ways for improvements. 

The paper is organized as follows: In Sec. II, the model we studied is described in detail. 
In Sec. Ill and IV, the simulation results are shown according to the removing strategies. 
In Sec. V, we give the conclusion. 



II. NETWORK TOPOLOGIES, ROBUSTNESS AND REMOVING STRATEGIES 



-V) model of random 
3j and the Newman- 



The network models of topology we studied are the Erdos-Renyi ( 
network the Barabasi- Albert (BA) model of scale-free network 
Watts (NW) model [a], which is a variant of the Watts-Strogatz (WS) model of small-world 
network [4]. For all the network models, we set vertex number N = 1024 and average degree 
(k) = 8. The edges are undirected and have no weight. Multiple edges are not allowed. In 
order to study the robustness, we ensure that the topologies are totally connected when they 
are intact. The ER model yields exponential degree distribution. Each vertex pair is linked 
by an edge with probability (k)/(N — 1). Since (k) > ln(iV), the ER networks are almost 
surely connected lj. The BA model features growth and preferential attachment. At any 
time step of a growth, the network from the previous step is connected. A new vertex is added 
with (k)/2 preferentially linked edges in the current step and the connectedness is preserved. 
If a BA network grows infinitely large, the probability of a randomly selected vertex having 
degree k is proportional to k~ 3 . The NW model starts with a regular structure, which is 
identical to the WS model. But unlike the WS model, the NW model builds shortcuts by 
randomly inserting pN edges, where p is the rate of the shortcuts. This way of building 
shortcuts eliminates the possible network fragment during the process of rewiring in the 
WS model As total connectedness is the special concern, we adopted the NW model 
instead of the WS model. With p changes from small to large, a NW network can undergo a 
transition from "large world" to "small world" . We choose p = 0.02 in the small-world phase 
without loss of generality. The pN shortcuts are inserted into a periodic square lattice with 
side length L = 32 and coordination number z — 8. As p is a small quantity, the average 
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degree is 8 + p ~ 8. Among the three models, the BA model produces scale- free degree 
distribution which is observed in many real networks The NW model is based on a 
square lattice, which is the topology that can reproduce some of the observed real Internet 
features 3(3]. The randomly inserted shortcuts are for the small- world property, though 
high clustering is absent. The ER model does not match real networks in nearly all aspects. 
Nevertheless, it has the significance that many of its properties can be obtained through 
probabilistic approaches. Note that the average path lengths of all the three models scale 
logarithmically with the number of vertices. It is interesting to compare the robustness of 
them. 

The robustness of the network models refers to the ability of maintaining total connected- 
ness when a portion of edges are removed. The total connectedness is defined as the effective 
communication between every pair of vertices, which is affected by the twofold effects of the 
edge removal. One is that the physical connection between a pair of vertices may be cut 
off. If the number of edges in a network is decimated to be less than N — 1, the network is 
certainly not totally connected. The other effect is related to the fact that a network usually 
performs some function, e.g., transmits some quantity. The edge removal may redistribute 
the quantity transmitting on the network and the vertices may lose communication due to 
congestion [ijj], e.g., the quantity from a source vertex is congested on some intermediate 
vertices and never reaches the target. The former effect is topological, while the latter is 
functional. Either of the two effects can destroy the total connectedness, and cause the 
breakdown of the networks. 

The topological breakdown of the networks occurs when the networks are just fragmented 
from one component into two during the edge removal, i.e., at least one of the vertices is 
isolated physically from the others 13] . This condition is more restrictive on the connectivity 
than the existence of the giant component, which was usually used in previous works {(], 3, 
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12j. We adopt this condition in favor of the emphasis on the completeness of a 



network and that the functional breakdown can be defined in a similar manner. 

We incorporate a simple traffic dynamics to define the functional breakdown of the net- 
works. At each time step, a certain amount p of quantity, which can be data, energy, etc., is 
transmitted between every pair of vertices, along the shortest paths. Transmitting along the 
shortest paths minimizes the delay of distributing the quantity. For the logarithmic average 
path length of the network models, the transmission scales well. The number of shortest 
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paths passing through a vertex j, which defines the vertex betweenness Bj of j [26|, deter- 
mines the load Lj of the transmitting quantity on j, Lj = pBj 14j. The load is handled by 
the vertex within its capacity, i.e., the amount of quantity that can pass through j is at most 
the capacity of the vertex. It would be efficient to set the capacity of j proportional to Bj, 
so as to meet the load. However, the calculation of betweenness needs global information 
and takes time O ({k)N 2 ) 23|, which is resource consuming. For simplicity, we set uniform 
capacity C for each vertex. As a result, the vertex with the largest betweenness -B max is 
the most likely to encounter overload. The largest generation rate p c of the transmitting 
quantity is C/S max . The rate p c marks the capacities of the networks. It was found that 
homogeneous networks have higher capacities than heterogeneous ones [l4 ] . As the purpose 
of the study is to compare the different aspects of the robustness rather than the absolute 
capacities, we normalize p with respect to p c for each individual network, i.e., the load of 
j changes to Lj = pBjC/B max , where p is the normalized generation rate of the quantity. 
The load Lj is initially less than C with < p < 1, but Lj could be larger than C as 
the vertex betweenness redistributes in the process of edge removal. If removing an edge 
leads to Lj > C for some vertex j, congestion occurs [141 ] . The excessive load accumulates 
continuously on j and the vertex is unable to communicate with other vertices. Similar to 
the topological breakdown, the total connectedness is destroyed. This situation defines the 
functional breakdown of the networks. One can see that the functional breakdown actually 
depends on the evolution of the ratio Bj/B max ; the study of the functional robustness is 
essentially a test of the response of vertex betweenness to edge perturbations. Though the 
dynamics is simple, some essential functional properties of the networks are reflected and 
we are able to study the relative importance of the functional robustness to the topological 
robustness under different values of the only parameter p. 

To study the two types of the robustness, we start from an undamaged network generated 
by one of the network models and remove the edges one at a time until the network breaks 
down for either the topological or the functional reasons. The robustness are studied in the 
contexts of different removing strategies. 

Two removing strategies are used: random failure and attack. The random failure strat- 
egy removes edges with uniform probability, which can be seen as a simple abstraction of 
the successive error in a communication network. The attack strategy removes edges in the 
descending order of their importance, which tries to model a sophisticated attacker who 
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knows the global information of a network and always targets the most important link. 
There are several quantities that can define the importance of an edge, such as edge be- 



tweenness and edge degree [10|, [22|. We choose edge betweenness as the definition of the 
importance, because betweenness directly measures the load of the transmitting quantity in 
our traffic dynamics. Analogous to vertex betweenness, edge betweenness is defined as the 



number of shortest paths passing through an edge 22] . Every time the most important edge 
is removed, we recalculate the betweenness to find the most important one in the remain- 
ing edges. If there are multiple edges having the same largest betweenness, we randomly 
select one for removal. Random failure and attack are usually studied together to obtain a 



6|. In the next two sections, we show the 



collective characterization of network robustness 
results according to the two strategies. 

III. THE RANDOM FAILURE STRATEGY 



In this section, edges are randomly removed from the networks. We first study the total 
breakdown, which is the combined effect of the topological breakdown and the functional 
breakdown. Figure [1] shows the cumulative distribution function CDF to tai of the total break- 
down as a function of the fraction / of removed edges for the network models. The results 
are shown for different values of p. 

When p = 0, there is no quantity transmitting on the networks. The total breakdown 
is in fact the topological breakdown, which is characterized by a critical fraction of 
removed edges. When / < /* opo , each of the networks is totally connected; when / > /* opo , 
each of the networks is split into two or more components, where one of the components 
could be the giant component. 



Thus 



For the ER model, almost every model realization is connected if (k) > In (AT) 
f*po=[(k)-ki(N)]/{k). 

For the BA model, a simple estimate of f c is sought under the framework of random 



graph with prescribed degree sequence [29j. Seeing that the BA model yields no assortative 



mixing 



28 



i.e., no degree correlation, we equate the BA networks to random scale- free 
networks. Although the BA model features network assembly and evolution instead of 
complete randomness, some qualitative result can be obtained by the estimate and confirmed 
through simulation. For a random network with arbitrary degree distribution P(k), the 
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FIG. 1: The cumulative distribution function CDFtotai of the total breakdown under the random 
failure strategy as a function of the fraction / of removed edges for (a) the ER model, (b) the 
BA model and (c) the NW model. When p = 0, the total breakdown is in fact the topological 
breakdown; when p > 0, the total breakdown is the combined effect of the topological breakdown 
and the functional breakdown. All the data is averaged over 10240 realizations. 



probability tt s of a randomly selected vertex being in a finite component of size s has been 
obtained in Ref. [271 ]. 



7T, 



(k) 



d s ~ 2 



dz 



(1) 



J 2=0 



where gi(z) = YX=o ( - fc+1 - > ^ ( - fc+1 - > z k and s > 1. Particularly, 7Ti = P(0). We now consider a 
scale-free random network, which has the same N, (k) and P(k) as the BA networks. With 
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the minimum degree larger than 2, gi(z) has no constant terms and tt s = 0. There are no 
finite components, all the vertices are connected in the giant component. It was reported 
that the giant component always exists for a random network with P(k) ~ k~ 3 when vertices 
are randomly removed 7j • The result is similar if edges are removed 12| . Thus the network is 
always split into the giant component and a finite component. Denote 7r s (/) the distribution 
of the size of the finite component after a fraction / of edges are removed. When / — > 0, the 
network is nearly not affected and 7r s (/) — > 0. Though 7r s (/) is small, various sizes of the 
finite components are probable. With the existence of the giant component, 7r s (/) ~ vr s (0) 
decays exponentially 27j]. We then neglect the higher order components and focus only 
on the components of size one and two. The new degree distribution and the new average 
degree after edge removal are P f (k) = E^ =fc P Oo) ($) (1 - f) k f ko ~ k and (k) f = (1 - f){k), 
respectively [7]. We calculate the ratio 



r(f) 



nif) (k) f Pf(o) 

(1-/) \^ =1 kP(k)f 



(2) 



f 2 (k)T,r =1 p(k)f* • 

This ratio is a monotonically decreasing function of /, which implies that in addition to 
the exponential decay of 7r s (/) when /—►(), the components of size two become even less 
probable than the components of size one as / becomes larger; in most cases, the network 
is split into the giant component and a finite component consisting of an isolated vertex. 
In the simulation of the BA model, we observed that in more than 99.8% realizations, the 
finite component is of size one. This fact allows us to estimate the critical fraction /* opo 
by the emergence of an isolated vertex, which has empty degree. If one such vertex can be 
sampled, the total connectedness is destroyed, i.e., fl° po satisfies P/ c to PO (0) 



A' ' 



For the NW model, the critical fraction f l ° v ° is obtained in a similar manner. Employing 
the expression of 7r s (/) [27|, we get 



fy-\k) 



d 



8-2 



(s-1)! 
(7s)!(A;)/ 6s+2 (l 



dz s ~ 2 
-fY 



(3) 



(s-l)!(6s + 2)! 

where the shortcuts in the model are neglected for convenience and the original degree 
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FIG. 2: The scaling relation between P / (0) and 1/N at f = f c & P° for the BA model and the NW 
model when edges are randomly removed. Symbols are the results of simulation, each of which is 
an average over 102400 networks. The two solid lines (They nearly overlap.) are the linear fits of 
the simulation results. The slopes of the lines are 0.991 and 0.980. 

distribution reads P(k ) = 5(k — 8). When s is large, 

KJ ' f] -5(1 -J)/ 6 - (4) 



7T,-l(/) 6 6 

This ratio tends to zero when / — > and increases monotonically until it reaches the 
maximum value 1 at / = f c — I, where f c is the critical percolation fraction [7|. As 
the total connectedness is more restrictive than the existence of the giant component, we 
claim that < fl° po < f c . In this region, vr s (/* opo ) drops exponentially, and again, each of 
the NW networks is split into the giant component and an isolated vertex, f* opo satisfies 
Py c to P o(0) = jr. Particularly, we get f l ° v ° ~ N~» for the NW model. The result is consistent 
that if N — > oo, then fj: opo — > and vr s (/* opo ) decays very quickly with increasing s. In Fig. [2, 
we plotted the scaling relation between Pjt opo (0) and 1/N obtained in the simulation for the 
BA model and the NW model. As the figure shows, the agreement between the theoretical 
estimate and the simulation is reasonable. The above calculations show that /*°p° is strongly 
correlated to the local connectivity of vertices. In the BA model, the majority of the vertices 
have degree four, while in the NW model, all the vertices have degree at least eight. Thus 
as Fig. [1] shows, the NW model is more topologically robust than the BA model and they 
both are more topologically robust than the ER model. 
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When p > 0, the functional breakdown can also happen. Intuitively, networks with 
larger p are more functionally vulnerable. We showed in Fig. [1] the cumulative distribution 
function of the total breakdown only for near capacity generation rates. When p = 0.95, 
the ER networks are very likely to break down when only a small fraction of edges are 
removed. For the same p, the NW model has slightly larger CDF tota i than the BA model, 
which suggests that while the NW model is more topologically robust than the BA model, it 
is less functionally robust. Note in the figure that while CDF to ta\ for large p increases fast in 
the region < / < 0.1 for all the network models, CDF tota i for p = increases much slower 
for the ER model or remains nearly zero for the BA model and the NW model. Thus in this 
region of /, the total breakdown is mainly determined by the functional breakdown. With p 
close to 1, the functional breakdown reflects the change in the network capacity during the 
edge removal. In Fig. [3], we plotted the average network capacity (p c ) / as a function of / for 
the network models. {{p c )f is normalized to the initial value (p c )o-) The capacity of all the 
network models decreases with the edge removal, because the shortest paths concentrate on 
the remaining edges and the betweenness of vertices becomes more heterogeneous. 

For the ER networks, the degree distribution remains exponential when edges are ran- 
domly removed, i.e., the ER networks are still ER networks after the edge removal. The 
vertex betweenness B of ER networks has a power-law relation B ~ k a with degree k 24| , 
hence 

C C 



(Pc); 



[fcmax(/j] 

C 



V R /f\ (l + a)[fcmax(/)]" 

2^3 D Ai) [fc max (/)] a+1 -l 



Ck mSuX (f) ^ 



where 5 max (/) is the maximum vertex betweenness and k max (f) is the maximum degree 



when a fraction / of edges are removed. The vertex betweenness a 
DN(N — 1), where D ~ ln(A)/ln((fc)y) is the average path length 



so satisfies J2j Bj(f) 
l|. Then we get 



Ck max (f)\n((k) f ) 
{Pc) f~ {l + a)N{N-l)hx{N) 1 } 



and the normalized average capacity 



(Pc)f 
(Pc)o 
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FIG. 3: The normalized average capacity (p c )// (Pc)o under the random failure strategy as a 
function of the fraction / of removed edges for the network models. The inset shows the linear 
relation between {p c ) f / (Pc)o and F(f) for the ER model, where F(f) = k ma _ x (f)ln[(k) (1 — /)]. 
All the data is averaged over 10240 realizations. The line in the inset is the linear fit. 

where k max (f) and / have an implicit relation [2| through the regularized gamma function 
P(a,x) that P [\_k max (f) + lj , (1 — f)(k)] = jj. The linear relation between (p c ) j / (p c ) 
and F(f) is also shown in Fig. [31 The average network capacity of the ER model decreases 
much more quickly than those of the other two network models, thus the ER model has the 
worst functional robustness. 

For the BA model, though there is a similar power-law relation B ~ k v in scale-free net- 



works 



25 



jScale-free degree distribution does not remain scale-free when edges are randomly 



removed [20||. Instead, we compare the trend of the the load distribution for the BA model 
and the NW model. The two network models both have heterogeneous load distributions. 
Hub vertices or shortcuts carry large load. The BA model also has heterogeneous connec- 
tivity while the NW model has homogeneous connectivity. When each edge is removed with 
equal probability, the hub vertices in the BA networks are more likely to be diminished than 
the non-hub vertices, while the shortcuts in the NW model are not biased. The BA net- 
works tend to have more homogeneous load distribution than the NW networks. As shown 
in Fig. [3], the average capacity of the BA networks decreases slightly slower than that of the 
NW networks. The BA networks are a little more functionally robust. 
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FIG. 4: The occurrence probability Ptopo of the topological breakdown under the random failure 
strategy as a function of the normalized generation rate p for the network models. The data is 
averaged over 10240 realizations. The dashed line is for Ptopo = 50%. The network model that 
corresponds to the smallest p at the point of intersection between Ptopo and this line has the best 
relative topological robustness. 

After studying the combined effect, we compare the dominance of the topological break- 
down and the functional breakdown. During the simulation, we measure the occurrence 
probability Ptopo of the topological breakdown. If this probability is larger than 50%, the 
topological breakdown is dominant, otherwise the functional breakdown is dominant. As a 
supplement to the absolute robustness studied above, the dominance is a sign of the rela- 
tive robustness, i.e., which breakdown is less likely to happen than the other. The result is 
shown in Fig. H] as a function of p. For all the network models, the topological breakdown is 
dominant when p is small, and the functional breakdown is dominant when p is close to 1. 
There is a shift in the dominance, which corresponds to a particular value of p. The smaller 
this value is, the better relative topological robustness a network model has. As the figure 
shows, the NW model has the best relative topological robustness as well as the best abso- 
lute topological robustness, and the BA model has the best relative functional robustness as 
well as the best absolute functional robustness. Though the ER model has both the worst 
absolute robustnesses, its relative robustness is intermediate. 
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IV. THE ATTACK STRATEGY 



In this section, edges are removed in the descending order of their betweennesses. For 
the sake of clarity, we first show the results for the ER model and the BA model, and then 
for the NW model. We study the same quantities as in the previous section. 

The cumulative distribution function CDF to ta\ of the total breakdown was plotted as a 
function of / in Fig. [5] for the ER model and the BA model. We first discuss the topological 
robustness, which corresponds to p = 0. The ER networks are very vulnerable to attack, the 
removal of only 0.5% edges almost surely destroys the total connectedness. The vulnerability 
is due to the lack of loops. With nearly no loops [l], the ER networks can be roughly seen as 
trees, i.e., there is only one self-avoiding path connecting a pair of vertices. The only path is 
at the same time the shortest path. Many shortest paths concentrate on the high betweenness 
edges, which are removed by the attacker with high priority. The loss of the edges prevents 
communication between vertex pairs easily, as there are no alternative paths. On the other 
hand, the BA networks are much more robust, they can afford the removal of nearly half of 

n 

the edges before almost surely break down. There is a strong correlation [10J in BA networks 
that (e) ~ k e , where (7b (e) is the betweenness of an edge e and k e is the product of the 
degrees of the two vertices connected by the edge. This correlation implies that at the early 
stage of the attack, when the characteristics of the networks are not affected too much, 
only the edges between vertices with large degree are targeted. The removal of these edges 
generally does not destroy the total connectedness and the BA networks can be reduced up 
to tree-like. The topological robustness of the ER model and the BA model has been studied 
by examining the size of the largest component in Ref. [10] . The authors focused on the 
effects of various attack strategies. Here we examine the total connectedness and provide 
explanations for the effects of the specific strategy. Moreover, the topological robustness is 
compared with the functional robustness. As shown in the figure, plots corresponding to 
empty load and high load almost overlap for both of the network models, indicating that 
the functional breakdown hardly occurs. We then investigate the average network capacity. 

The normalized average capacity (p c )/ / (Pc)q was plotted as a function of / in Fig. E]for 
the ER model and the BA model. When the edge with the highest betweenness is attacked, 
the shortest paths that originally pass through the edge take detours. The highest edge 
betweenness is dispersed and the network capacity is increased. This is the case for both 
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FIG. 5: The cumulative distribution function CDF tota i of the total breakdown under the attack 
strategy as a function of the fraction / of removed edges for (a) the ER model and (b) the BA 
model. All the data is averaged over 10240 realizations. The plots corresponding to empty load 
and high load almost overlap for both of the network models. 

of the network models. The capacity of the ER networks monotonically increases during 
the whole process, indicating that the networks are free of the functional breakdown. The 
capacity of the BA networks keeps increasing until a certain value of /, after which the 
shortest paths revert to collect on the remaining edges and the capacity is decreased. As 
shown in the figure, the capacity is boosted up to nearly 8 times as the initial value at 
/ = 0.35 and then drops. We observed in the simulation that only in few realizations 
does the capacity drop below the initial value, thus the BA networks are almost free of 
the functional breakdown. Note that the growth and decay of the capacity was studied in 



Ref. |18J with the purpose of finding the optimum transmission efficiency. At each time 
step, the author chose to remove the edge with the largest weight, which is the product 
of the betweennesses of the two vertices an edge connects, and intentionally avoided the 
disintegration of network. Though there are some variations in the models, our result is 
qualitatively no different. 
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FIG. 6: The normalized average capacity (p c ) f I {Pc)o under the attack strategy as a function of 
the fraction / of removed edges for (a) the ER model and (b) the BA model. All the data is 
averaged over 10240 realizations. 



An interesting point worthy of note is that the attack strategy does not fulfill its name for 
the BA model. By comparing Fig. [5](b) with Fig. [T](b), we find that the BA model is more 
robust under the attack strategy than under the random failure strategy for any value of p. 
This result coincides with the moderated vulnerability of networks in terms of edge [l^] and 
is in contrast to the higher vulnerability of the BA model under vertex attack than error 

n 

when the existence of the giant component is used as the criterion of breakdown |6J . 

The cumulative distribution function CDF tota i of the total breakdown and the normalized 
average capacity (p c )/ / (Pc) are shown as functions of / for the NW model in Fig. [7J The 
figure can be interpreted as follows: Starting from a complete NW network, the attacker first 
targets the shortcuts, because they collect the most shortest paths. After all the shortcuts 
are removed, the rest of the network is a two dimensional L x L periodic lattice with each 
vertex connected to its neighbors and next neighbors. The structure becomes completely 
homogeneous and the network capacity is boosted. Thus we see the first peak of the network 
capacity at / = fx — 2p/(k) = 0.005 in Fig. [7J(b). The attacker continues by randomly 
targeting an edge. We show a diagram in Fig. [HJ^a) to help the interpretation. The lattice is 
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put on a x-y plane and Fig. E^a) shows a slice along the x direction. For every edge in the 
figure, there are L — 1 parallel counterparts in the y direction. The L edges form a group 
that if any edge in the group is removed, the others are removed in the next few steps as 
well, because shortest paths that originally pass through the whole group concentrate on the 
remaining members, raising the betweennesses of them. Suppose that edge e% in the figure 
is removed, the network capacity drops until all the counterparts of the edge are removed. 
The heterogeneity of the raised betweenness is alleviated and we see the second peak of the 
network capacity, which is L edges away from the first peak. However, the capacity is not 
fully restored, as edge e% and their counterparts collect larger number of shortest paths 
than before. Thus two peaks follow, each at an interval of L edges. With the 3L edges gone, 
the lattice is no longer periodic and the center area in the x direction collects the highest 
betweenness. Then the following attack removes another 3L edges in this area and the 
network is split. For the most of the model realizations, six peaks of the network capacity 
are observed, which conforms to the arguments above. There is a special case that allows 
the observation of one more peak. As depicted in Fig. Mjo), it is possible that a shortcut 
is placed between vertex a and b. This kind of shortcut spans less lattice distance than the 
abundant next-neighbor links, thus collects less shortest paths and is removed very late in 
the attack. The existence of such shortcut brings higher betweenness to a group of L edges 
in the y direction, e.g., all the shortest paths passing through a and c collect on the edge e§ 
in the figure. Removing this group of edges induces the seventh peak. The seven peaks are 
pointed out in Fig. [71(b) by arrows. We see in the figure that the network capacities on the 
peaks are larger than the initial value, thus only when the capacities in the valleys drop below 
the initial value does the network functionally break down, i.e., CDF tota i shown in Fig. UJa) 
resembles stairs for p ^ 0. When p = 0.95, nearly 90% model realizations functionally 
break down before / = 0.005. The situation is only a little better than the ER model and 
much worse than the BA model. For the topological case p = 0, there is a plateau in the 
figure which corresponds to the emergence of the seventh peak. We estimate the emergence 
probability as (^^jP sc {l — P sc ) pN ^ 1 ~ 7.3%, where P sc = is the probability of the 
occurrence of the special shortcut and we only consider the first order of this probability. 
The result is supported by Fig. [7] that the topological breakdown takes place just before the 
seventh peak of the network capacity, at / = /V — fx + 6L/((k)N/2) = 0.0519, for about 
92.7% model realizations. The topological breakdown for the rest of the model realizations 



17 




FIG. 7: The simulation results of the attack strategy for the NW model: (a) The cumulative 
distribution function CDFf t a i of the total breakdown as a function of the fraction / of removed 
edges for different values of the normalized generation rate p; (b) The normalized average capacity 
(Pc) f I (Pc)q as a function of /. All the data is averaged over 10240 realizations. The arrows point 
out the seven peaks of the network capacity. The dashed lines are guides to the eye: the topological 
breakdown (p = 0) takes place just before the seventh peak for about 92.7% model realizations. 

takes place after L more edges are removed. The NW model has intermediate topological 
robustness between the ER model and the BA model. 

The relative robustness of the network models for the attack strategy is studied in the 
same way as in the previous section for the random failure strategy. The result is shown 
in Fig. [91 The ER model has the worst absolute topological robustness, and also the worst 
relative topological robustness, as for any value of p, no functional breakdown takes place. 
While the BA model has the best absolute topological robustness, the relative topological 
robustness is still bad, nearly the same as the ER model. For the NW model, we know that 
the topological breakdown only takes place when at least 6L edges are removed, but the 
functional breakdown could happen if network capacity drops below the initial value. The 
figure shows that for p > 0.2, the topological breakdown hardly takes place. Thus the NW 
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FIG. 8: A NW network after all the shortcuts are removed, which is a two dimensional periodic 
lattice with coordination number 8, is put on a x-y plane. Both the arcs and lines represent edges, 
and the circles represent vertices, (a) A slice of the lattice in the x direction is shown. Every edge 
in the complete lattice has equal betweenness. If e\ is removed, the shortest paths going from the 
left-hand side to the right-hand side have fewer choices of edge than before. Thus ei and es collect 
larger betweenness than the others and are the immediate targets in the following attack. Note that 
ei, e2 and e% each represent a group of L edges in the y direction, (b) The edge is such a special 
shortcut that it spans less lattice distance than the next-neighbor links and collects fewer shortest 
paths. Though the betweenness of the shortcut is low, the betweennesses of the neighboring edges 
are increased, e.g., the shortest paths connecting a, c and a, d both pass through 6, the edges e§ 
and e§ are biased. Different from e$, the edge es is in the y direction. In addition to the six peaks 
in the x direction, one more peak of the network capacity could be seen. In fact the effect of the 
shortcut e4 is weak that the counterpart of es which connects to a in the y direction is not removed 
before the lattice is split. For clarity, some next-neighbor links are omitted. 

model has the best relative topological robustness. 
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FIG. 9: The occurrence probability Pt opo of the topological breakdown under the attack strategy as 
a function of the normalized generation rate p for the network models. The data is averaged over 
10240 realizations. The dashed line is for Ptopo = 50%. The network model that corresponds to the 
smallest p at the point of intersection between Pt opo and this line has the best relative topological 
robustness. 

V. CONCLUSION 

In this paper, we have studied the topological robustness and the functional robustness of 
several network models under two strategies of edge removal, using the total connectedness 
as the measure. For each removing strategy, we have examined the combined effect and the 
relative importance of the two types of robustness. Through the study of the combined effect, 
we have found out the network topology which is the most robust in a specific environment, 
e.g., the NW model is the most robust in the environment of the random edge removal, while 
the BA model is the most robust in the environment of the edge attack. Through the study 
of the relative importance, we have known with evidence how to strengthen a network effi- 
ciently, e.g., as the NW model has the best relative topological robustness, improvements on 
the functional robustness should be emphasized, such as increasing the capacity of vertices; 
the BA model has bad relative topological robustness under the edge attack, thus enhance- 
ments on the connectivity such as building redundant links are appropriate. These results 
can have applications in designing and optimizing artificial networks, such as implementing 
a robust P2P network where a connection between peers has a constant probability to fail. 
There are also some extensions for further studies. The traffic dynamics in our model is far 



20 



from realistic, more elements abstracted from real-world traffic can be incorporated. More- 
over, networks are not limited to carrying traffic. It is an open question how the interaction 
between the topologies and different dynamics influences the network robustness. 
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